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Quantum mechanical calculations based on an extended Heisenberg model are compared with 
ferromagnetic resonance (FMR) experiments on prototype trilayer systems Ni7/Cu n /Co2/Cu(001) 
in order to determine and separate for the first time quantitatively the sources of the temperature 
dependence of interlayer exchange coupling. Magnon excitations are responsible for about 75% of 
the reduction of the coupling strength from zero to room temperature. The remaining 25% are due 
to temperature effects in the effective quantum well and the spacer/magnet interfaces. 
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The coupling of two magnetic layers through a non- 
magnetic spacer layer (interlayer exchange coupling, 
IEC) has received considerable attention recently. The 
physics at T=0 can be described quite generally in the 
framework of a quantum interference model P, e.g. 
the dependence on the spacer thickness. Within this the- 
ory, the period, the amplitude, and the phase of the cou- 
pling are well understood today. However, most of the 
experiments on magnetic multilayers and their IEC have 
been undertaken at ambient temperature, not only for 
easier feasibility but also in view of technological appli- 
cations. In this Letter we evidently show in theory and 
experiment that the major effect is due to magnon excita- 
tions, whereas band structure effects are less important. 

Much experimental and theoretical work focused on 
the oscillating dependence of the IEC on the spacer thick- 
ness resulting in ferromagnetic (FM) and antiferromag- 
netic (AFM) coupling of the magnetic layers [4J. As for 
the temperature dependence there is convincing exper- 
imental evidence on several IEC systems that the cou- 
pling strength decreases with T from its T=0-value to 
almost zero at T — Tq in metallic systems Jl, @, 0| and 
that it shows an effective T 3 / 2 -dependence [5(. However 
relatively little theoretical work has been carried out con- 
cerning the temperature dependence of the IEC over the 
full temperature range. The detailed ab-initio investi- 
gations which focused on Fermi edge softening effects 
are restricted to low temperatures by the assumption 
that the spin-resolved band structure does not change 
significantly [g, |£|. Working at all temperatures, Mills 
and coworkers investigated the influence of spin waves 
and predicted the experimentally found T 3 / 2 -dependcncc 

In this Letter we present a novel approach which com- 
bines a quantum mechanical treatment of magnetic mul- 
tilayers at all temperatures with experiments on the pro- 
totype IEC systems Ni7/Cu„/Co2/Cu(001) (subscript: 
number of layers). Up to now a quantitative separation 
of the possible sources of the T-dependence (spacer ef- 



fects like softening of the Fermi surface or reduction of 
the spin asymmetry of the reflection coefficients vs. mag- 
netic layer effect, i.e. excitation of magnons) has not been 
achieved The key benefit of the current approach is 
the possibility to "switch on and off" magnon excitations 
in our model. Thus, by performing one-parameter fits of 
the microscopic interlayer exchange coupling with and 
without magnons, we can separate the two main sources 
of the temperature dependence of the IEC. We find that 
the magnetic contribution exceeds the spacer part by a 
factor of 3 in the range from w 0.15 • Tq to To 

The theoretical description of T-dependent IEC exper- 
iments often relies on the classical, macroscopic Landau- 
Lifshitz (LL) equations. This continuum's approach is 
usually used together with an expansion of the free en- 



ergy to obtain the IEC energy by defining [111. \Y1\ 



2Jintcr(T) = Ftec(T) = Tjj(T) - F n {T) 



(1) 



Although the LL approach allows for an examination of 
the overall T-dependence of the IEC, it cannot distin- 
guish between the different contributions, i.e., it cannot 
yield © in the form F IEC (T) = F spacer (T) + T magn .(T). 

Figure Q] shows FMR experimental evidence (evalu- 
ated with LL) that the magnetic layer effect gives in- 
deed a significant contribution. According to the the- 
oretical works by Bruno [l[ and Edwards et al. 0] the 
T-dependence should be more pronounced for a higher 
number of spacer layers while it should be independent 
of the spacer thickness within the spacer theory proposed 
in Ref. [9( . The observed stronger decrease of the 4-ML- 
spacer (AFM coupled) system as compared to the 5-ML- 
spacer (FM coupled) system suggests that there are other 
than just spacer mechanisms at work. Moreover one sees 
a large difference between the measured Ji n tcr-curve of 
the 4-ML-spacer film and the theoretical curve according 
to the spacer theory suggested in Ref. 

Ferromagnetic resonance is a powerful tool to investi- 
gate the IEC experimentally [HJ ■ This technique probes 
the uniform spin wave mode w(q = 0) of a magnetic 
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FIG. 1: Temperature dependence of Jinter normal- 
ized to Jintef (T=0) for the trilayer IEC system 
Ni7 / Cu n /C02 /Cu(OOl) for three different spacer thick- 
nesses ncu 1J]. The data for 9 ML (full squares) are taken 
from Ref. The dashed line corresponds solely to a spacer 
effect contribution for the nc u = 4 ML sample. 
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FIG. 2: Resonance field for Ni 7 /Cu(00f) and 
Cug/Co2/Cu(001) as a function of the orientation of the 
external magnetic field (circles). The angle 9b is measured 
with respect to the film normal. Lines are theoretical fits, 
Ni: 5=1, K 2 = 3.0 /iBkG, go = 4.5 /i B kG, J = 30 meV, Co: 
5=2.5, K 2 = -20.25 /z B kG, g = 1.95 /i B kG, J = 4.1 meV. 



sample. An external field Bo is tuned for a given probe 
frequency i\f = uj(q = 0)/2n until resonance occurs at 
B les (T,9B ), with 9b being the angle between the mag- 
netic field and the normal to the film plane. This res- 
onance field B rcs (T,9B ) at which uniform (q = 0)-spin 
wave modes with the energy E sw — hv^i of the prob- 
ing microwaves are excited in the magnetic system is 
the crucial quantity to connect experiment and theory. 
With the method we propose it is possible to compute 
the spin wave energy £J sw (q = 0,Bo,T) and read off 
the resonance field given the (experimental) parameter 
i>hi, which was set to 9 GHz in our calculations 
For further details about the experimental procedure we 
refer to Ref. 12| and give here a short outline of the 
experimental procedure. After the preparation of the 
Cu n /Co 2 /Cu(001) single film already capped with n Cu 
spacer layers, the angular dependence of the resonance 
field B^° s (T,6bo) was determined at two temperatures 
(55 K and 294 K). Then the temperature dependence of 
the resonance field B?°{T, Bq ) with T E [55 K, «300 K] 
at 0b o = 90° was measured. The same was carried out 
after the deposition of the top nickel layer for the cou- 
pled system Ni7/Cug/Co2/Cu(001) and for comparison 
also for the single film system Ni7/Cu(00T). 

The Hamiltonian of the effective two-layer system we 
use to describe IEC reads 



H — — JgSigSja — ffq^BBpS 

<ij>a ia 

— ^ JjSfaSjff — 7^ K-ZaS. 
<ij>a{3 ia 



(2) 
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The first term describes nearest-neighbor exchange cou- 



pling J a between spin moments Si a and S Jct at sites i 
and j within the same layer a (=1,2). The film thick- 
ness is implicitly included in the parameters J a which 
are chosen such that the magnetic moments of the two 
monolayers described in (2) equal that of the consid- 
ered multilayer films at room temperature, respectively 



Ni 7 



410 K, T, 



Co 2 



400 K). The second term con- 



tains an external magnetic field Bo in arbitrary direc- 
tion with the spectroscopic splitting factors g a (taken as 
2.2 for Ni and Co) and the Bohr magneton /.iq. The 
spacer layer is not considered explicitly. The IEC me- 
diated by the spacer is described by the parameter Jj 
(third term) . The fourth and fifth term represent second 
order lattice anisotropy and dipolar interaction, the lat- 
ter leading to shape anisotropy. Ki a and go a are micro- 
scopic anisotropy parameters, Si az is the z-component 
of Si a and is perpendicular to the film plane, r\j is 
the distance between sites i and j, and dipolar coupling 
across the spacer is neglected. The shape anisotropy fa- 
vors in-plane orientation of the magnetization, the lattice 
anisotropy in-plane (K2a < 0) or out-of-plane (K2a > 0) 
orientation. Our Hamiltonian is similar to that used in 
Ref. [1 51 ] for the investigation of magnetic reorientation 
and anisotropy. We also point out that the local-moment 
aspect of inherently itinerant magnets like Ni and Co has 
been demonstrated before: we mention here Heisenb erg - 
like spin wave dispersions of an itinerant band model [16j , 
a high degree of the d-band moment localization (> 90% 
for < T < Tq) in an LDA + many-body approach 
to bulk Ni [l7|, and the successful employment of the 
Heisenberg model with ab-initio exchange parameters for 
the calculation of the Curie temperature of Ni [18( . 

We summarize shortly how we solved ([2]) for the single- 
magnon Green function (GF) ((Sj~ a ; Sj~p)) and refer to 
Refs. [IH, [3, [2(| for details. In order to treat the 
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FIG. 3: Resonance field for the coupled system 
Ni 7 /Cu9/Co 2 /Cu(001). All parameters for the Ni and 
Co subsystems are the same as in Fig. [2] 
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FIG. 4: Resonance field at 9b = 90° as a function of tem- 
perature for Ni7/Cu9/Co2/Cu(001). Same parameters for Ni 
and Co as in Figs. [2] and [3] Dashed line: the T-dependence 
due to the softening of the Cu Fermi surface [l[ was included. 



anisotropy terms properly we rotate the coordinate sys- 
tem self-consistently such that the new z-axis is parallel 
to the direction of the magnetization. In the subsequent 
equation of motion approach it is important to include 
all four combinations of the GF {{S^~; S~h~)) as was 
first suggested in Ref. [2(j. This ensures correct soften- 
ing properties of the uniform spin wave mode, namely its 
vanishing at the reorientation field for Ki > and up to 
a critical field strength for K 2 < 0, respectively. Higher 
order GFs stemming from the Heisenberg exchange and 
from the dipolar interaction are decoupled using the 
Tyablikov (RPA) method. In the case of the dipolc- 
dipole interaction only the uniform (q — > 0) contribution 
is considered as the nonuniform terms are negligible com- 
pared to contributions from the much larger Heisenberg 
exchange. For the local iiVterms an RPA decoupling 
fails. It was shown in Ref. [19( that the Anderson-Callen 
decoupling yields quantitative agreement with QMC re- 
sults for the field induced magnetic reorientation transi- 
tion. This is also the case for our improved theory which 
combines the proposal of Ref. with the multiple GF 
approach by Pini et al. [2o| . Different from the theoret- 
ical treatment based on the Landau-Lifshitz equations, 
now with the quantum Heisenberg model a self-consistent 
evaluation of the magnetization and the spin wave exci- 
tation spectrum E sw (q, B , T) is carried out. 

In our theory the anisotropies Ki a and g$ a influence 
the system solely via the effective anisotropy field given 
by the temperature-dependent term 



K 2a {T) = (S za )(T) (2K 2a C a (T) - Dg 0a ) 
(S(S + l)-(S 2 za )(T)) 



C a [T) = 1 



2S2 



(3) 
(4) 



Here the T-independcnt quantity D is some number de- 
pending on the lattice geometry. Note that our effec- 
tive anisotropy field K 2a corresponds to the quantity 
M e g commonly used within a Landau-Lifshitz descrip- 
tion of FMR experiments [12j. Furthermore we exploit 



(S 2 za )(T) = S(S + I) - {S za )(T)[l + 2<p a (T)} with the 
average magnon occupation number <p a [13]. The T- 
dependence of K 2a goes beyond a mere proportionality 
to (S za )(T) due to the expectation value (Sl a )(T). 

Figure [2] shows the comparison between the B ics (8b )- 
curves from theory and experimental data for a 
Cu 9 /Co 2 /Cu(001) and a Ni 7 /Cu(001) film system at 
two different temperatures. In both cases the effective 
anisotropy favors the magnetization to lie within the film 
plane. There is quite good agreement at both tempera- 
tures over the whole range of angles 9b for both films. 

It is important to note that for a fixed temperature 
the B tcs (9b ) data alone do not allow one to determine 
the anisotropy parameters Ki a and g§ a independently of 
each other [see Eq. Q]. But taking into account both 
the -Bres(^B ) data and B rcs (T) at fixed #b o =90° the pa- 
rameters can be separated because of the temperature 
dependent factor C a (T) in Eq. ([3]). Indeed it is still 
possible to accurately fit the experimental results with 
one set of (T-independent) parameters {K2 a , go a ) for Ni 
and Co, respectively, over the whole temperature range. 
Furthermore in both cases the values of ga a S lie slightly 
above the result of an explicit evaluation of this quantity 
assuming point-like dipoles on the lattice sites for the 
given geometry [12] (goS = 3.81/XBkG). The conclusion 
we draw is that the T-dependence of the effective mag- 
netic anisotropy is solely due to spin wave excitations 
which manifest themselves in the T-dependence of Q 
rather than due to T-dependent Ki a ^ go a which would 
effectively describe nonmagnonic sources as thermal ex- 
pansion or phononic interactions. This means that we 
only have one fit parameter, J/, for the IEC system. 

Figure [3] shows the resonance fields for such a film 
system. The two branches at a given T correspond to 
the acoustic (lower field branch) and optical (high field 
branch) excitation modes, respectively. The optical mode 
lying above the acoustic indicates antiferromagnetic cou- 
pling between Ni and Co. The missing experimental data 
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FIG. 5: Temperature dependence of interlayer exchange cou- 
pling in Ni7/Cu9/Co2/Cu(001). Symbols stem from fits of the 
experimental resonance field at 9 = 90° at the given temper- 
atures. Open circles: T-dependent calculation. Filled circles: 
calculation at T = 0. Lines are guides to the eye. 



merits on prototype Ni7/Cu ra /Co2/Cu(001) interlayer ex- 
change coupled films. Employing a microscopic theory 
based on an extended Heisenberg Hamiltonian we stud- 
ied the temperature dependence of the IEC by fitting 
experimental data on the FMR resonance field in two 
steps, first with T-dependent calculations and then with 
T = 0; thus, "switching off" magnon excitations in the 
system. The result is that magnon excitations are most 
important (ss 75%) for the temperature dependence of 
the interlayer exchange coupling. Band structure effects 
(e.g. smearing of Fermi surface) are present but of minor 
importance. As a result we note that a theoretical treat- 
ment of IEC at T = is of limited use when interpreting 
the IEC of real multilayers at temperatures relevant for 
technological applications. 

Financial support by the DFG (Sib 290) is gratefully ac- 
knowledged. 



for the optical mode at T = 55 K are due to a low inten- 
sity (oscillator strength) of the resonance signal [12J . All 
parameters except for Jj stem from the fits of the single 
(uncoupled) Ni and Co films in Fig. [2] 

In addition we carried out fits at a fixed angle 6 — 90° 
for a variety of temperatures, the results of which are 
shown in Fig. 2J The variation of J/ with T represents 
the spacer contributions. From comparison with the re- 
sults obtained by additionally taking Fermi surface soft- 
ening into account we conclude that sources other than 
that dominate the spacer contribution, presumably the 
T-dependence of the reflection coefficients. 

In order to separate quantitatively the spacer and the 
magnetic T-dependence of the IEC we fitted in a sec- 
ond step the same -B ros -data but setting the magnon oc- 
cupation number in our model artificially to zero, i.e., 
we "switched off" the magnons. The whole temperature 
dependence of the IEC including magnon excitations is 
then accounted for by the fit parameter jj . From the 
comparison between both resulting curves of Ji{T) and 
^ji^ ^ n pig Pg] one can reaa ; ff the weight of each con- 
tribution to the T-dcpcndcncc of the IEC. One concludes 
that the magnonic part (ss 75%) exceeds the spacer part 
(w 25%) significantly, however the latter is not negligible. 

As for the variation of both contributions with the 
number of magnetic and non-magnetic layers one can 
state the following general trend: the fewer spacer lay- 
ers and the thinner the magnetic part (i.e. the smaller 
Tc), the more important spin waves become for IEC. Of 
course, for bulk-like systems of Fe or Co with T^ ulk > 
1000 K the magnon part will be (possibly much) less sig- 
nificant at room temperature or less. However, for Tq 
below « 700 K, as for thin nano-film systems or Ni films 
or manganite multilayers 2l| , the spin wave contribution 
to the IEC temperature dependence becomes dominant. 

In conclusion we proposed a novel approach that com- 
bines a quantum theoretical treatment of the magnetic 
properties of multilayer film systems with FMR experi- 
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